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ABSTRACT 

We present molecular dynamics simulations of the motion of a single rod in a two- 
dimensional random static array of disks. For long rods the mean-squared displacement 
of the center-of-mass shows a cage effect similar to that observed in supercooled liq- 
uids or dense colloidal systems. We have determined the time-dependence of the non- 
Gaussian parameter for different rod lengths. It is found that the long-time regime is 
strongly non-Gaussian even at length scales of the order of 10-15 times the rod length, 
thus showing the heterogeneity of the dynamics at such length scales. 



§1 INTRODUCTION 

The Lorentz gas is a simple model in which a particle diffuses in a random 
configuration of static obstacles. Initially introduced by Lorentz as a model for 
the electrical conductivity in metals (Lorentz 1905), it can be used as a rough ap- 
proximation for the dynamics of a light atom in a disordered environment which 
presents a much slower dynamics. It has been shown that, for high densities of 
obstacles, it qualitatively shows many of the features of the dynamics of super- 
cooled liquids or dense colloidal systems, such as a transition from an ergodic 
phase of non-zero diffusivity to a non-ergodic phase with no diffusivity (Bruin 
1972, Alder and Alley 1978, Gotze et al. 1981a,b, Masters and Keyes 1982). De- 
spite its simplicity, this problem is highly non-trivial even in the simplest case 
where the diffusing particle and the obstacles are modelled as hard-spheres. In 
particular, diffusion constants and correlation functions are non-analytical func- 
tions of the density of obstacles (Bruin 1972, Gotze et al. 1981a, b, Binder and 
Frenkel 1990). In this work we present an investigation by means of molecular 
dynamics simulations on a generalization of the Lorentz gas: the diffusion of a 
rod in a random configuration of static obstacles. This picture can be under- 
stood as a cartoon of the dynamics of a linear molecule in a porous medium or in 
a colloidal suspension of heavy particles. Though simulations have been done in 
two dimensions, we do not expect qualitatively different features in the general 
three-dimensional case. 

The paper is organized as follows: In Section 2 we present the model and 
give some details of the simulation. General dynamical features are shown in 
Section 3. A study on the non-Gaussianity of the long-time regime is presented 
in Section 4. Conclusions are given in Section 5. 



§2 MODEL AND DETAILS OF THE SIMULATION 

The rod consists of N beads of mass m and radius a equal to that of the 
obstacles. The rod length is therefore 2Na. The density of obstacles is defined 
as p = Aobs/^box' with -Wobs the number of obstacles and -Lbox the length of the 
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square simulation box. Beads and obstacles interact via a soft-sphere potential 
V(r) = e(<r/r) 12 . The latter is truncated and shifted at a cutoff distance of 2.5a. 
In the following, space and time will be measured in the reduced units a and 
(cr 2 m/e) 1 / 2 respectively. 

A set of 100 different realizations of the rod was propagated in a given con- 
figuration of the obstacles. Periodic boundary conditions were imposed. The 
set was equilibrated at temperature T = e/fce by coupling it to a stochastic 
heat bath. After the equilibration, a run of 10 s time steps, with a step size of 
0.01, was performed at constant energy. In order to improve the statistics, this 
procedure was repeated for 6-10 different configurations of the obstacles and the 
obtained results were averaged over the corresponding runs. 

§3. TRAJECTORIES AND DYNAMICS 

In the following, all the results will be presented for a fixed density of ob- 
stacles p = 6 x 10~ 3 . An extensive investigation for a wide range of N and p 
will be presented elsewhere. Figure 1 shows typical trajectories for two char- 
acteristic rod lengths N = 5 and 50. They respectively correspond to rods 
of dimensions comparable and much longer than the average distance between 
nearest-neighbour obstacles, d n . n . « 7. For N = 5, the trajectory presents a 
structure of blobs and jagged channels. Thus, the rods are able to penetrate 
often in regions of high local density of obstacles, where they are trapped for 
some time, forming the blobs, while the jagged channels correspond to paths in 
regions of lower local density, where the collisions are less frequent. 

A very different polygonal-like structure is obtained for N = 50, consisting 
of nearly longitudinal motions in the tubes formed by the obstacles, connected 
by vertexes corresponding to regions where the rod leaves the tube and enters 
into a new one. 

Figure 2a shows the mean-squared displacement of the center-of-mass of 
the rod (r 2 (t)) divided by the time t. Brackets denote ensemble average and 
r(t) = |R(i) — R(0)|, with R(t) the position of the center-of-mass of the rod at 
time t. At short times the rod does not feel the presence of the obstacles and 
(r 2 (t)) shows a quadratic time-dependence, as expected for a ballistic motion. 
For small N a sharp transition to the long-time regime (r 2 (t)) ~ t is observed, 
while for long rods a crossover regime, which extends over 1-2 decades, develops 
at intermediate times. This corresponds to the well-known cage effect observed 
in supercooled liquids (Kob and Andersen 1995, Sciortino et al. 1996, Caprion 
and Schober 2000, Mossa el al. 2000, Colmenero et al. 2002) or dense colloidal 
systems (Puertas et al. 2002, van Megen 2002, Weeks and Weitz 2002): the 
particle is temporarily trapped by the cage formed by their neighbours and as 
a consequence, the collisions within the cage result in a time-dependence of the 
mean-squared displacement that is slower than the one in the ballistic regime. 
Finally, the particle escapes from the cage and begins to diffuse. It is worthwhile 
to remark that in the present case the breaking of the cage is not achieved by 
cooperative motions between neighbouring particles, as in the case of liquids or 
colloids, since the configuration of obstacles is static. Instead, the trapping of 
the rod is of purely entropic nature. It must also be noticed that the observed 
cage effect is much less pronounced than in supercooled liquids or dense colloids. 
Higher densities of obstacles or longer rods would be required for observing a 
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stronger effect. 



§4. NON-GAUSSIANITY 

In Einstein's random walk model, particles move under the effect of colli- 
sions with the others. Between two consecutive collisions, the particle moves 
along a straight line in a random direction that is independent of its previous 
history. If the time and distance between two consecutive collisions are much 
smaller respectively than the observational time and lenght scales, the van Hove 
self-correlation function takes a Gaussian form (see for example Hansen and Mc- 
Donald 1986). This result is exact for an ideal gas and for an harmonic crystal. 
It is also valid in the limit t — ► 0, where atoms behave as free particles. At 
time and length scales comparable to those of the dynamics of the collisions, 
i.e., at observational scales sensitive to heterogeneities of the system, deviations 
from Gaussianity are expected. They can be quantified by the second-order 
non-Gaussian parameter C(2{t), which for a two-dimensional motion is given by 
&2(t) = [(r 4 (t)) /2(r 2 (t)) 2 } — 1. It is straightforward to see that for a Gaussian 
function a2(t) = 0. 

Figure 2b shows the TV-dependence of the non-Gaussian parameter. While 
for short rods a broad structure is obtained, for long ones a peak develops, 
which grows and shifts to longer times for increasing N. A qualitatively similar 
behaviour has been reported for other very different systems as a function of 
different control parameters, as temperature for supercooled liquids (Kob and 
Andersen 1995, Sciortino et al. 1996, Caprion and Shober 2000, Mossa et al. 
2000) or volume fraction for colloid-polymer systems (Puertas et al. 2002). 
However, in contrast to these systems, we find that the peak does not grow in a 
completely monotonic way, and that for short rods (N = 5, 10 in figure 2b), it 
seems to be located at longer times that for longer rods. 

From a comparison with figure 2a, it is deduced that for long rods (N > 10) 
the time t* for the maximum of 0:2 (i) is located at the time interval connecting 
the end of the cage regime and the beginning of the long-time linear regime, 
i.e., the breaking of the cage leads to a high probability of jumps much longer 
that the size of the cage, resulting in a strongly non-Gaussian -or heterogeneous- 
dynamics at that time scale. Similar results have been reported for other systems 
(Kob and Andersen 1995, Sciortino et al. 1996, Caprion and Shober 2000, Mossa 
et al. 2000, Colmenero et al. 2002). The size of the cage r c can approximately 
be estimated as r c = (r 2 (t*)) 1 / 2 . We find a value of r c w 35, independent of 
N within the error bar. This value is notably larger than the typical distance 
between neighbouring obstacles d n n « 7. However, it must be noted that the 
latter perhaps is not the most appropiate length for characterizing the caging 
dynamics of the rod trapped by the neighbouring obstacles. In contrast to a 
liquid-like equilibrium configuration, which shows an approximately homogeneus 
structure, a random configuration is characterized by clusters of close particles 
and big holes. For the density p = 6 x 10~ 3 , the size of the holes, rh, is typically 
between 40 and 70, which introduces another characteristic length. The size 
of the cage should lie at some value between d n . n . and rh, consistent with the 
estimated value of r c . 

In the case of short rods (N < 10), where no intermediate cage regime 
is observed, the maximum of ct2(t) does not seem to be located at the sharp 
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transition between the ballistic and the long-time regime, but at a later time. 
The reason for this difference with the large-TV case might be that, as argued 
in Section 3, for these lengths the rod has a high probability of exploring the 
holes, and at the same time is able to penetrate in regions of high local density, 
leading to non-Gaussian dynamics at length scales of the order of the hole size. 

The most striking feature in figure 2b is that in the time limit of the simula- 
tion, a2(t) has not decayed to zero but remains finite, especially for long rods, 
i.e., the long-time regime is still non-Gaussian at the length scale covered by the 
simulation. This effect can be studied by investigating the q-dependence of the 
incoherent intermediate scattering function F s (q,t) = (exp(— iq- [R(t) — R(0)])), 
where brackets denote ensemble and angular average. 

Figure 3 shows F s (q, t) for N = 50 and several values of the wavelength 
A = 2w/q corresponding to length scales much larger than the size of the rod. 
A two step-decay separating the cage and long-time regimes is not observable 
for this range of A, as in the case of supercooled liquids, where the plateau in 
F s (q,t) is only observable at q- values not too far from the position of the first 
maximum of the structure factor. 

In Gaussian approximation in two dimensions, the incoherent intermediate 
scattering function is given by F s (q,t) = exp(— (r 2 (t))q 2 /A). Therefore, Gaus- 
sianity in the long-time regime is only fullfillcd if the time dependence of the 
long-time decay of F s (q,t) can be well described by an exponential function 
exp(— t/r(q)), since at that time scale (r 2 (t)) ~ t. Figure 3a shows the cor- 
responding fits to exponential functions. It must be stressed that, usually, an 
amplitude factor A{q) < 1 must be included in the fit to take into account the 
strength of the long-time regime. While in other cases this is an important pa- 
rameter, for the present case setting A(q) = 1 did not change the results from 
the fit within the error bar, since the value of F s (q,t) is very close to unity 
at t ~ 10 3 , the time for the beginning of the long-time regime (see figure 2a). 
From figure 3a it is clear that only for wavelengths A beyond « 2000 -i.e., about 
20 times the rod length- the decay of the curve can be reasonably described in 
terms of an exponential function, and even for such large scales non-Gaussianity 
cannot be excluded, since the final decay is not accessed within the time window 
of the simulation. 

A better description can be achieved (figure 3b) by fitting to a Kohlrausch- 
Williams- Watts (KWW) function, exp[— (i/r(g))' 3 ], a phenomenological equa- 
tion often used in the analysis of long-time relaxations of complex systems 
(Phillips 1996). It must be noted that these type of fits are not fully satis- 
factory: as can be seen in figure 3b, the final part of the long-time decay relaxes 
slower than a KWW function. In any case, the obtained values of (3, different 
from unity (see inset in figure 3b), again confirm the non-Gaussian nature of the 
long-time linear regime at large length scales. 

§5. CONCLUSIONS 

We have carried out molecular dynamics simulations for a simple model of 
the motion of a rod in a random static environment. The dynamics of long 
rods show similar features to those observed in supercooled liquids and dense 
colloids, namely a cage effect at intermediate times between the ballistic and 
long-time regimes, due to the trapping of the rod inside the tube formed by the 
neighbouring obstacles. 
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The time-dependence of the non-Gaussian parameter has been calculated 
for different rod lengths. It is found that it does not decay to zero within the 
time window of the simulation, which for long rods cover length scales of the 
order of 10-15 times the rod length. Therefore, the dynamics is heterogeneous 
on such observational scales. The origin of the observed non-Gaussianity might 
be the above mentioned presence of big holes in the random configuration of 
obstacles, that would lead to a finite probability of jumps much longer than the 
average, corresponding to processes where collisions provide a nearly longitudinal 
instantaneous velocity to the rod and at the same time the rod finds a long free 
"corridor" between the obstacles. An investigation of the dynamics of the rod 
in an homogeneus liquid-like configuration of the obstacles, where big holes are 
not present, could shed new light on this question. Work in this direction is in 
progress. 
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Figure 2: (a): Mean-squared displacement of the center-of-mass divided by t for 
rod lengths N = 5, 10, 15, 20, 30, 40, and 50 (top to bottom). Points mark the 
time t* for the maximum of the non-Gaussian parameter. 

(b): Time-dependence of the non-Gaussian parameter for different rod lengths. 




Figure 3: Incoherent intermediate scattering function F s (q,t) for N = 50 (solid 
lines) and different values of A = 2n/q. Dashed lines are fits to an exponential 
(a) and to a KWW function (b) . The inset in (b) shows the dependence on A of 
the stretching exponent (3. 
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